Relationship between spiral and ferromagnetic states in the Hubbard model 

in the thermodynamic limit 
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We explore how the spiral spin(SP) state, a spin singlet known to accompany fully-polarized 
ferromagnetic (F) states in the Hubbard model, is related with the F state in the thermodynamic 
limit using the density matrix renormalization group and exact diagonalization. We first obtain an 
indication that when the F state is the ground state the SP state is also eligible as the ground state in 
that limit. We then follow the general argument by Koma and Tasaki [J. Stat. Phys. 76, 745 (1994)] 
to find that: (i) The SP state possesses a kind of order parameter, (ii) Although the SP state does 
not break the SU(2) symmetry in finite systems, it does so in the thermodynamic limit by making 
a linear combination with other states that are degenerate in that limit. We also calculate the one- 
particle spectral function and dynamical spin and charge susceptibilities for various ID finite-size 
lattices. We find that the excitation spectrum of the SP state and the F state is almost identical. 
Our present results suggest that the SP and the F states are equivalent in the thermodynamic limit. 
These properties may be exploited to determine the magnetic phase diagram from finite-size studies. 

PACS numbers: 71.10.-w,71.10.Fd 
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I. INTRODUCTION 

The itinerant ferromagnetism has been one of the most 
fundamental problems in condensed matter physics, and 
we are still some way from a full understanding of the 
problem. For the ferromagnetism in the Hubbard model, 
a simplest possible model for interacting itinerant elec- 
trons, there are various studies, which |are pioneered by 
Gutzwilleiu, Hubbarda, and KanamoriB. For the Hub- 
bard model on finite-size systems, it recently extended 
to rigorous proofs that spin-independent Coulomb in- 
teraction can indeed result in fully-polarized ferromag- 
netic ground states without degeneracy -for appropriate 
conditions and/or appropriate modelsaO. Apart from 
these, there is a body of numerical studies for finite sys- 
tems, where various authors have shown that finite Hub- 
bard models with appropriate band fillings and bound- 
ary conditions have fully-polarized ferromagnetic ground 
state&El. 

On the other hand, for infinite systems, the itiner- 
ant nature of electrons makes the problem interesting 
and subtle. The most important reason is the exis- 
tence of the spiral-spin(SP) state, a spin-singlet state 
that has the spjii-apin correlation length as large as 
the sample sizell3"ll3, which is known to accompany the 
fully-polarized ferromagnetic (F) state in various finite- 
size itinerant models. By accompanying we mean SP and 
F states are close in energy in finite systems, where the 
ground state changes between them when the boundary 
condition is twisted. 

The SP state does not seem to be an accident of the 
Hubbard model, since SP states are found not only in 
electron systems with strong charge-charge interactions, 
but also in those witfi spin-dependent risteractions such 
as double-exchangec3 or Kondo latticecll models. So we 



believe that it is a general feature of strongly correlated 
systems for the F state to be accompanied by the SP 
state. It is rather surprising, since the latter has the to- 
tal 5* = while the former has S = S'max- Thus it is 
desirable to resolve this question to understand the fer- 
romagnetism in itinerant systems in the thermodynamic 
limit. Thus the problem we have to clarify is the relation 
between the SP and F states in the thermodynamic limit 
(do those states degenerate or coalesce, etc), which we 
address in this paper. 

First, we show that the SP state is eligible as the 
ground state (i.e., degenerate with the F state in en- 
ergy in the thermodynamic limit) when the F state is 
the ground state. We check this for one-dimensional (ID) 
t-t' Hubbard model with infinitely large coulomb interac- 
tion taken as a typical example and calculate the energies 
of the SP and F states as a function of the inverse sys- 
tem size with the density matrix renormalization group 
(DMRG) method. 

Next, we discuss whether the symmetry of the SP state 
is broken in the thermodynarnie limit, following the ar- 
guments by Koma and Tasaki23 on symmetry breaking 
and finite-size effects in quantum many-body systems on 
a lattice. We show that (i) the SP state possesses a kind 
of order parameter, (ii) although the state does not break 
the global SU(2) symmetry in the finite system, it does 
so in the thermodynamic limit by making a linear com- 
bination of other states that are degenerate in that limit. 

We have studied excitation spectrum of the SP state 
by calculating one-particle spectral function, dynamical 
spin and charge susceptibilities for the SP state in vari- 
ous finite-size models, including ID t-t' Hubbard model, 
Tasaki's fiat band model, and Hubbard ladder. We have 
found that the excitation spectrum of these two states in 
the finite size system are almost identical. 
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The present results suggest that we should regard the 
SP state as identical with the F state although they have 
the total spins in the opposite extremes {S = for the SP 
state, S = iSmax for F) in the finite size system. These 
properties can be exploited to determine the magnetic 
phase diagram from the calculation of the total spin of 
the finite-size system in an unambiguous fashion. Oth- 
erwise we have to worry about not only the total spin, 
but the spin-spin correlation to identify SP states in de- 
termining the phase diagram. 

The paper is organized as follows. In Sec. II, we first 
confirm that the SP and F states have degenerate energy 
for L —f oo{L: sample size) in that the energy per site is 
identical between them. DMRG is used to calculate the 
difference in the total energy, AE = Esp — Ep, between 
the SP state and the F state as a function of system size 
L for the t-t' Hubbard model with U = oo as a typical 
example. The extrapolate AE(L — > oo) is always finite 
(or zero depending on the boundary condition even in 
the thermodynamic hmit), i.e., AE per site does indeed 
vanish. It implies that when the F state is the ground 
state, the SP state is also eligible as the ground state if 
these states are translation invariant. Then we move on 
to the discussion on symmetry breaking of the SP state 
along the Koma-Tasaki's argument. In Sec. Ill, we cal- 
culate the single-particle spectral function and dynamical 
spin and charge susceptibilities for various models. Our 
results suggest that the SP state and the F state are 
identical state in the thermodynamic limit. 

In Sec. IV, we take an appropriate combination of the 
boundary condition and the number of electrons to make 
the SP state higher in energy than the SF state, and 
the phase diagram for the one-dimensional t-t' Hubbard 
model, Tasaki model, and 2-leg Hubbard ladder model 
is obtained with exact diagonalization of small systems. 
We have found that the phase boundaries rapidly con- 
verge to those for larger systems obtained with DMRG. 
We may expect that it can be exploited to the determi- 
nation of the phase diagram of 2D or higher dimensional 
systems. A summary of the present study is given in Sec. 
VI. 



II. DISCUSSION ON THE THERMODYNAMIC 
LIMIT 

A. The energy difference between the SP state and 
the F state in the thermodynamic limit 

ID t-t' Hubbard model is one of the simplest models 
that is thought to exhibit ferromagnetism for sufficiently 
strong electron-electron interaction. Miiller-Hartmann 
suggested that in the low-density limit the, ground state 
of i-i'|_Hjibbard model is ferromagneticcj. Daul and 
NoackOa have carried out highly accurate DMRG cal- 
culation to conclude that there is an extensive ferromag- 
netic phase in the phase diagram. Since the t-t' Hubbard 



model is numerically tractable, we study the energy dif- 
ference of the SP and F states by taking this model with 
L/ = oo as an example. 

The t-t' Hubbard Hamiltonian is given by 

L 

n = -tJ2J2'^clc^+l,. +il.C.) 
2—1 a- 

L L 
1=1 (T i=l 

where cf^. creates an electron at site i with spin cr(=t 
,1), t is the nearest-neighbor hopping, and t' the next 
nearest-neighbor hopping, U the Hubbard repulsion, and 
Uia- = cl^Cia. In this section we take t = 1, t' = 0.2 and 
U = oo. 

As we can check by a simple exact diagonalization, the 
ground state of the system has always S'tot = for any 
value of U at least for system sizes up to 12 sites when a 
periodic boundary condition(PBC) is adopted. This is to 
be contrasted with the phase diagram obtained by Daul 
and Noack with an open boundary condition(OBC), in 
which a wide ferromagnetic region is found for t' = 0.2 
and large U . 

Let us first confirm that the SP and F states are degen- 
erate in the sense that their difference in the total energy 
behaves as 

AE = E^singict - £^F ^ finite, 
i.e., the energy per site vanishes, 
AEjL 0, 

where i^singict is the lowest energy within the S'tot = 
sector, namely the energy of the SP state. 

We take OBC, because the DMRG becomes most ^CjfUi- 
rate for this condition. As Daul and Noack showedHEJ, 
the ground state of t-t' Hubbard model is ferromagnetic 
for t' — 0.2 and [/ — > oo in OBC at least up to 50 sites. To 
estimate AE, we must calculate the energy of the ground 
state of the S'tot = sector which has higher energy than 
the ferromagnetic ground state. Therefore, we add the 
term AStot to the original Hamiltonian (Q) to selectively 
shift the states with higher total spin to higher energies 
by turning on A > 0, while conserving the SU(2) symme- 
try. Here we set A = 1. In Fig. n^, we show the spin corre- 
lation function, (S'f S|) [i = L/2, j = L/2 + 1, • • • , L - 3) 
of the spin singlet state for L = 32 sites, n — 0.75 as 
an example. We can see that the spin correlation wave 
length is as large as the system size. Namely, the spin 
singlet state is the SP state. 

Since DMRG is a variational procedure, the energy of 
the SP state is an upper bound, while the energy of the F 
state is calculated exactly because the F state does not 
feel the on-site coulomb interaction. Hence we overes- 
timate AE. To minimize this overestimation, we must 
calculate the energy of the SP state as accurately as pos- 
sible. This is the reason why we have set \J = oo and 
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t' = 0.2. Namely, the exclusion of double occupancies 
reduces the Hilbert space drastically, while a small value 
of \t'\ reduces the truncation error. 

Using the finite-size algorithm in DMRG, we calcu- 
lated up to L = 48 for the density of electrons n = 0.5, 
sweeping the system about 20 times to improve the wave 
function. We store the density matrix at each step to 
construct good initial vector for each super-block diag- 
onalization. We have kept up to m = 500 states per 
block at each step, where the convergence is checked by 
comparing the results for m — 300,400,500. The trun- 
cation error is smaller than 10~^, which is small enough 
to enable us to extrapolate AE{L — » oo). 

In Fig. ^(solid line), we show the results for AE = 
Esp — Ep as a function of inverse system size l/L for 
electron density n = 0.5. We can see that all the points 
for L > 24 fall upon a linear dependence on l/L, from 
which we can extract AE{L — > oo). For all the densities 
studied, AE{L — > oo) indeed remains finite 10~^). 

In order to check whether the result is not an accident 
for U — oo, we can introduce the effect of large but finite 
U as an (antiferromagnetic) exchange interaction J > 0, 

J^Pgi [si-Si+i-in.n.+i'jPG, (1) 

added to the U ~ oo t-t' Hubbard model. Here Si = 
\ "^13-1 ^\i3'^ Pi^ii spin operator, where a is the Pauli 

matrices and Pq denotes the Gutzwiller projection oper- 
ator. In Fig. ^ the results for J ~ ±0.2 are superposed. 
We can see that AE does not change drastically even if 
we introduce the effect of finite U . 

This kind of AE may depend on the boundary condi- 
tion, which can be indeed the case with finite systems. 
For instance, Kusakabe and Aokill3 have shown for a fi- 
nite two-dimensional Hubbard model that Nagaoka's fer- 
romagnetic state (one hole in the half-filled band with 
U ^ oo) alternate with an SP state as the boundary 
condition is changed from periodic to anti-periodic. In 
other words a level crossing takes place between the two 
states as an Aharonov-Bohm magnetic flux is introduced 
to twist the boundary condition. Such a boundary effect 
can affect the 'energy gap' even in the thermodynamic 
limit in general. 

Thus we compare in Fig. ^ the ordinary OBC (the 
same as Fig. ^, solid line) and the boundary condition 
in which we further turn off the nearest-neighbor trans- 
fers [t) at either end of the system in the t-t' model for 
AE{L — > oo) at n = 0.5. We can see that AE{L oo) 
becomes zero (negligibly small) in the latter case. There- 
fore, we can not exclude a possibility for which the SP 
state has a lower energy than the F state in some appro- 
priate boundary condition. In other words, the problem 
which state has a lower energy is a very subtle problem. 



B. The definition of the ground state 

As we have seen, it is meaningless in the thermody- 
namic limit to identify the ground state by studying small 
'finite' differences in the total energy between the candi- 
dates, so that a totally different point of view is required. 
The mathematical definition of the ground state in the 
thermodynamic limit is reviewed in Koma and Tasaki's 
articlec3 and here we follow them. For simplicity, we fo- 
cus on the case of PBC. 

We first recapitulate the definition. Let A be an arbi- 
trary local operator that acts on a finite number of sites. 
We define p as 

p^^^ = Ji'?.^A(A) 

for each A, and we call p as 'state' in the thermodynamic 
limit. Here A is a finite {Lx Lx ■ ■ ■ L x L) d-dimensional 
hypercubic lattice while Z"^ is an infinite d-dimensional 
lattice, and we have defined 

Pa(---) = TrHA[(- • Opa], 

where the trace is taken over the Hilbert space Ha 
spanned over A, and pA is an arbitrary density matrix 
on Ha. We also assume a Hamiltonian Ha that is local 
as defined by 

Ha = ^ hx ■ 

where is a local component such as = 

J2cr{cl:a'^x+i,a + il-C.) + Unx^nxi for the Hubbard model. 
We then define the ground-state energy density Eq as 

£0= lim inf 7^($a|Ha|«'a), 
AT^" *AeHA, 11*11=1 |AP 

where $ is a state in H\ and |A| is the norm of A. A 
state uj in the infinite- volume limit is said to be a ground 
state if it satisfies 

U!{hx) = Eo 

for arbitrary x. 

In the present context, we have seen in the previous 
section that the numerical results for the Hubbard model 
imply that the energy per site of the SP state and the 
F state is identical in the thermodynamic limit for any 
boundary condition, especially PBC. In this boundary 
condition the system has a translational invariance, so 
that an SP state in a finite system should have 

r(x)|SP)A =exp(z7)|SP)A 

under the translation T{x) which translates the state by 
X. We can expect that the invariance holds in the ther- 
modynamic limit as well, so that 
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(SP|/i,„+,|SP) = {SP\T{x)h,,T-\x)\SP) 
= (SP|/ix„|SP) 

= £o, 

and we can conclude that the SP state in PBC is also the 
ground state in the thermodynamic limit. 

C. Symmetry breaking in the spiral spin state 



On the other hand the order parameters, when not 
squared, satisfy 



,(SP|0i''|SP)A =A (SP|0f |SP)a = 0. 



Thus we have 



" ||(01)^-^|sp)a| 



One important question is: does the SP state, despite 
its being spin singlet, accompany a symmetry breaking 
in the thermodynamic limit? For simplicity, we consider 
the ID case. From the numerical results for the spin- 
spin correlation in finite systems, a natural quantity that 
is expected to become nonzero (in PBC) is 



-y 



cos(0, -0,)(S, •S,)sP >0, 



(2) 



where 9j = 2'Kj/L. We can then introduce an order pa- 
rameter, 



o 



(1) 



-sot 



ieA 



o 



(2) _ 



y ^ - sin e^Sf + cos 9iS', 



ieA 



(3) 

(4) 
(5) 

(6) 
(7) 

(8) 



where = S"" ± iSV. Three operators 0^^\ o'^^\ o'^^\ 
which obey a Lie algebra, satisfy the assumption required 
for Theorem 2.4 in Ref. |2^. The expectation values of the 
squared o'j^^ is 



,(SP|0^^^^|SP 



2 
3l2 



cos( 

ijSA 



9j){Si ■ Sj)sp 



1 



sin(0,-0,)(5+5-)sp. 



We can always take |SP)a to be real in PBC since 
TYa is real. Sf can also be represented as a real ma- 
trix, so that {S^S~)sp is a real number, and we have 
(5/S'j")sp — {Sj'S~)sp so that the last line above van- 
ishes. From the inequality eq.(E|) we have 



A(SP|Oi'''|SP)A = A(SP|OriSP)A 

2 



-,(2)2, 



3L2 ^ 

ijGA 



cos(0,-0,)(S, •S,)sP >0. 



as the 'low-lying state' in the Theorem 2.4 in Ref. |2^, 
which asserts that a symmetry breaking can occur when 
some low-lying excitations whose energies approach that 
of the state in question like 1/L mix with it. Here "ifj^ 
does have an energy that approaches to Esp as 

/\,f2 

KM/f ,i/A*f ) -Esp\< const. X ( — ). 



We also have, according to Theorem 2. 

lim lim(S^,O^^^S^) = const. > 

k-*ooA1Z ^ 



where, 



"A = 



1 



V2k+1 



A j 



M=l 



Here, let us introduce the translation operator T(x), 
which satisfies T{x)cl.^T^^{x) = c].^^^^. Since 

Tix)0+T-\x) - 1 5] exp(-z0.jT(a:)5+T-i(x) 

xaeA 

exp(-i6i^JS'++^^ 

XaeA 

T{x)'i>t^ = exp(iAf ^a;) exp{ij)'i'j^ , where 7 is a constant. 



Thus, for o 



+ 



,A/-1\ 



Therefore, we have 



lim lim 

k^ooATZ \U\ 



Al"a; 



xen 



appearmg m eq. (J^j, 

A(vI/f|e-'^^+-5++,Jvi/M-i)^ 

AK^|e-''==r(x)o+T-i(a;)K-i)A 
g(-ie^+iAfe^-i(M-i)e^ +17-17) 

XA«|0+K-^>A 

a«|<K-^)a. 



' ' xen 



= const. > 0, 

where is a finite subspace of A whose size is . Hence, 
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lim lim^ V(Si|^+|S'= 



A/ 



const. > 



Namely, although the SP state itself does not break the 
symmetry, it does so in in the thermodynamic limit by 
letting {^^f} mix with it in that limit. In other words, 
an arbitrarily large but finite region in an SP state has a 
broken SU(2) with a finite magnetization. 

If we use identities S~ = gA^i^' [^^ j "S" ] — Oi 
(SPIS"- = 0, 



=A (SP|(Oa 
= 0. 



M) 



-M 



Thus both of the SP state and {-^f} that mixes with it 
in the thermodynamic limit turn out to be orthogonal to 
the F state for finite systems. 

Finally, let us discuss whether S is the ground 
state in the thermodynamic limit by checking whether 
(Sj/ia-jS) = Sq holds for arbitrary x. Since 



.(SP|(o;^)^^^/i.„+.(o+)*^|SP)A 

=A {SP\{0]^rTix)K,T-\x){0+r\SP)j 



'|SP>y 



, hx"^^) = £o does indeed hold for arbitrary x in the 
limit of L — > oo. This means S is the ground state. 



III. THE EXCITATION SPECTRUM OF THE 
SPIRAL SPIN STATE 



Having shed light in the context of Koma-Tasaki's ar- 
gument, we may now conjecture that the SP state and the 
F state are identical in the thermodynamic limit. Thus 
let us have a closer look at the SP state as compared with 
the F state in three one-dimensional models, i.e., the t-t' 
Hubbard model, Tasaki's flat band model, and the two- 
leg Hubbard ladder. To determine whether the SP and 
F states are identical in the thermodynamic limit, we 
must show that the relation a;sp(^) = ujsf(A) holds for 
arbitrary local operator A. Although it is impossible to 
confirm this for all the possible A, if the F state and the 
SP state are identical in that limit, the excitation spec- 
trum of these state in finite size system are expected to 
be almost identical. To confirm this, we look specifically 
at the single-particle spectral function and the dynami- 
cal spin and charge correlation functions of the SP state. 
Indeed we will show that since the nature of the SP state 
in finite systems is almost identical to that of the F state, 
we cannot distinguish the SP state from the F state by 
looking at these quantities for sufficiently large system. 



A. t-t' Hubbard model 

In the t-t' Hubbard model where the electron transfer 
extends to next nearest neighbors (t'), the ground state is 
Stot = for any value of U in PBC at least for sizes up to 
12 sites as mentioned in the previous section. Let us first 
look at the spin-spin correlation function, {^c\SiSj\^G)^ 
to confirm that the S'tot = ground state for large enough 
U > Uc is the SP state. We show in Fig. |j the re- 
sult for 10 electrons in a 12-site ring (two holes) with 
[/ = 20 or 40. For U = 20, a short-range antiferromag- 
netic spin-spin correlation is observed for {SiSi+i) < 0. 
By contrast, for U = AO (or larger, which falls upon the 
ferromagnetic region in the phase diagram; see Fig. |l^ 
below), the spin-spin correlation has a wave length as 
large as the system size, which is the same behavior of 
the lowest 5tot = state in OBC (Fig. |). 

A more detailed nature the states is encapsulated in 
the single-particle spectral function given by 



1. 



1 



LU±{Eo-n)- iO 



itl^G), (9) 



where A'^{A ) denote the electron addition (removal) 



spectrum with ^'^^ 



(7t)^ I*g) andi;o the 



ground state and its energy, respectively. We have nu- 
merically calculated this quantity. In Fig. ^, we show 
the result for the parameter values adopted in Fig.^ 
with U = 40. We can see that the Luttinger relation, 
kp — 7rn/2 (with n = 10/12 = 0.83 here), for the spin un- 
polarized electrons is violated in favor oikp = nn, which 
would be expected for fully spin-polarized fermions. This 
behavior in the single particle spectral function is con- 
sistent with the nearly ferromagnetic nature of the SP 
state. Such a behavier was first found for the t-t'- J model 
by Eder ancL-OhtaEll, and then given an interpretation 
subsequentlyES. 

The /c-dependence of the spectral function may be 
analyzed hv evoking the argument by Doucot and 
Wen(DW)ca, who studied the infinite-[/ Hubbard model 
on a two-dimensional square lattice with two holes to find 
a trial state that gives an energy lower by 1 /L^ than Na- 
gaoka's ferromagnetic state. The key idea of DW is based 
on the following intuition. Holes behave as free fermions 
for on-site interactions when the background spin state is 
ferromagnetic (or nearly so). Assume that fully polarized 
electrons take an open-shell configuration (in which there 
is a degeneracy in the free-electron configurations) in k- 
space in PBC. If the background spin state is changed 
from F to the SP spin texture, a fictitious gauge field cor- 
responding to half flux quantum is generated. This shifts 
the /c-points by half the fc-point spacing, which lowers the 
kinetic energy because the polarized electrons now take 
a closed-shell configuration. 

It is a nontrivial question whether the DW's trial state 
may be applied to reproduce the spinless-fermion-like be- 
havior in A{k,uj) found here for the ID system. Let us 
assume, for simplicity, that the holes hop in a rigid spin 
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background, although in the original paper DW consid- 
ered the holes dressed by spin waves. The energy of the 
hole in the t-t' Hubbard model is then given as 



-2t cos {k±^ + 2t' cos 2 (^fc ± 



(10) 



where k{= 2ttN/L,0 < N < L,iV is an integer) is 
wave number and ± corresponds to the sign of the fic- 
titious flux. We can see that the fitting in Fig. ^ 
is remarkably good. To be precise, the transfer inte- 
grals in the original DW are taken to be effective ones 
(t — > t cos{tt / L) , t' — > t' cos{2Tr / L)) to take care of the 
reduction in the transfer between slightly twisted spins, 
but the present result is better fitted with these reduction 
factors omitted (dashed curves). 

In the addition spectrum an almost dispersionless band 
of low intensity peaks is seen (at around E = 2 in Fig. 

If the ground state were the fully polarized F state, 
then such a band is expected, since we can add an oppo- 
site spin at any fc-point. A remnant of such a band again 
suggests the nearly ferromagnetic nature of the SP state. 

A next question is whether the SP state persists for 
more than two holes. In Fig. ^, we show the spin-spin 
correlation for 6 electrons on a 12-site lattice (six holes, 
or n = 0.5; quarter filled) in PBC for U = 4,6. We can 
see that for U = 6, the spin-spin correlation is again spi- 
ral. To investigate whether DWs picture is valid for such 
a high hole doping, we have calculated the single-particle 
spectrum for J7 = 6 in Fig. 0. The spectrum for the SP 
state is fitted well by the energy dispersion defined by the 
eq. (|l^) , again t and t' not reduced, even though the hole 
concentration is as large as rih — 0.5. This is surprising, 
since the assumption that holes are nearly free would be 
valid only for small enough doping. 

We next question whether the SP state is connected 
adiabatically to the antiferromagnetic state as we de- 
crease the Hubbard U. We show the ground-state energy 
as a function of 1/t/ for 6 electrons on a 12-site ring in 
Fig. ^. We can clearly identify a level crossing appearing 
as a cusp around [/ ~ 5, which indicates that a transition 
(rather than an adiabatic connection) occurs within the 
S = space from the antiferromagnetic phase to the SP 
phase. 

We now turn to the dynamical spin and charge corre- 
lation functions 



N{k,Lo) = ilni($G|Affc— — — 

TT U} + [Lq — 7t) — iO 



5(fc,c.)=^ilm(<i>G|5," 



1 



uj + {Eo-n)- iO 



where N^. and S" are Fourier transform of (n^ — n) and 
S°'{a = X, y, z), respectively. In Fig. ^, we show the re- 
sults for the F state in APBC and the SP state in PBC 
for 10 electrons in 12 sites with [/ = 40. We can see that 
the behavior of these functions for the SP state is sur- 
prisingly similar to those of the F state even for L = 12. 



We can expect that they will become identical in the 
thermodynamic limit. 



B. One-dimensional Tasaki model 

TasakJ^ proposed a Hubbard model on a special class of 
lattice structures for which the lowest energy band is dis- 
persionless (flat). He proved rigorously that the ground 
state is ferromagnetic for arbitrary interaction strength 
{0 < U < oo) when the flat band is half-filled. 

To clarify whether metallic ferromagnetism can be re- 
alized in models haviag a non-half- filled, nearly-flat band, 
Sakamoto and Kuball3 investigated with DMRG the Hub- 
bard model on a chain of triangles, which may be thought 
of as a realization of Tasaki's model in ID. They obtained 
results which suggest that the system exhibits metallic, 
fully-polarized ferromagnetism for sufficiently strong ip-. 
teraction. On the other hand, Watanabe and Miyashitata 
found that an SP state is the ground state at least for sys- 
tem sizes up to L = 12 for appropriate conditions, i.e., 
sufficiently large U , n > 0.5, even number of electrons 
and PBC. 

The Hamiltonian is given by 



L/2 



i—1 (7 
L/2 

+ 4+l,a)c2»,ff + H.C.] 

i—1 cr 

L/2 L 
i—1 (7 i—1 

where t is the transfer for the bottom bond in a triangle, 
at for the bonds connecting bottom and top sites and (3t 
the site energy of a top site. Since a unit cell contains two 
sites, the single-electron spectrum consists of two bands 
with the dispersion relation, up to a constant, 



2 cos A: ± a/ 4 cos^ k + (Sa^ _ 4/3) cos fc + /^^ 



8a^ 



The lower band is then flat when /3 = — 2 is satisfied. 
Hereafter in this subsection, we take i = l,a = 2, /3 = 
(for which the lower band is dispersive) and L even. 

While Watanabe and Miyashita have shown that SP 
states appear for more than quarter-filled bands (0.5 < 
n < 1.0), we find here that the SP state is the ground 
state of S'tot — sector for sufficiently large U in the 
region n < 0.5 as well at least for L < 12. In Fig. |l^, 
we show the spin correlation for 6 electrons on a 12-site 
lattice in PBC for C/ = 4 and U = 6. We can see that 
for the case oi U — 6 the spin correlation length is in- 
deed as large as the system size. We have also calcu- 
lated the single-particle spectral function, eq.(H), where 



we take 77^. 



If we combine 
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them, they should contain the contributions from the two 
bands. Here is a Fourier transform of 7^!^. In Fig. 
[n] , we show the result for 6 electrons on 12 site lattice, 
where the spectrum is fitted by a two-band extension of 
eq.®, 



e± = tcos(fc± y) ±tWcos2(fc± y) + 8cos(/c± y) + 8. 

Fig. |ll| shows that the DW picture (almost free elec- 
trons hopping in a twisted spin background) is surpris- 
ingly good even for the two-band case. 



C. Two- leg Hubbard ladder 

Ferromagnetic ground state in ID has also been ob- 
tained for ladder systems, where two chains are connected 
with an inter-chain transfer, t^. Liang and PangEj sug- 
gested that for ti^/t = 1 and U — 00, the F state is one of 
the ground states for n/j < 0.22, wHle the ground state 
is spin singlet for rih > 0.4. Kohnai3 presented rigorous 
results for C/ = 00 2-leg Hubbard ladder in the limit of 
large inter-chain hopping {t±/t 00). He also studied 
the case for finite values of t±/t with DMRG to obtain 
the phase diagram. His results are consistent with that 
of Liang and Pang. These studies assumed OBC. 

On the other hand, it is intriguing to study whether 
an SP state appears in the ladders as well for PBC, and 
if so, whether A{k,u!) exhibits a DW-like behavior. The 
Hamiltonian is given by 



W = -iEE(4-C.+ l,a<.+H.C.) 

i—1 a, a 

L L 
-i-L^^(4l<TCi,2<T +H.C.) + U ^^n^a^fliai, 



i—1 a 



n—1 a 



where a(— 1, 2) labels the two legs of the ladder. We set 
t = 1, f7 = 00. 

We have calculated the intra-chain spin-spin correla- 
tion {Sf^Sj^) (Fig.^, and found that the ground state is 
indeed SP. We have also calculated the spectral function, 

c^2cr (creating an 



J 

--klo 



^k2 
„t 
-■kla 



(an- 



eq. (|^) , where we now take 7^ 

electron in the bonding band) or 7^^^ — u^j^^ — 1.^,2^. 
tibonding). In Fig.|l3|, we show the result for two holes 
(14 electrons on a 8 x 2 ladder). We fit the spectrum by 
the dispersion. 



-2t cos 



±t. 



(12) 



Again the picture of two holes hopping in a twisted spin 
background gives an accurate description. 



IV. MAGNETIC PHASE DIAGRAM FROM 
FINITE-SIZE STUDIES 



Our results in the previous sections suggest that we 
should regard an SP state as being ferromagnetic rather 
than non-magnetic even though ^tot = 0. So we fi- 
(ll)ially come to the problem of how to determine the mag- 
netic phase diagram from finite-size studies in the light 
of the SP state. As we have stressed, the ground state 
of finite systems in a certain boundary condition always 
has Stot = no matter how U is strong in some one- 
dimensional models as exemplified by the t-t' Hubbard 
model and the Hubbard ladder. Even when a ferromag- 
netic state appears for finite systems, the magnetic region 
in the phase diagram can shrink—fpr some system size 
as shown by Sakamoto and Kubc£3 for the ID Tasaki's 
model. This is due to the appearance of the SP state, 
which are encountered if we assume PBC for even num- 
ber of electrons, or more generally, if we assume a bound- 
ary condition for which the fully-polarized electrons take 
an open-shell configuration in the ground state. 

Since the SP state may be regarded as ferromagnetic 
as elaborated in previous sections, we have then to dis- 
tinguish the SP state from nonmagnetic states by cal- 
culating the spin correlation, etc. Here we propose that 
this difficulty can be readily overcome by taking an ap- 
propriate combination of the boundary condition and the 
electron number for which the SP state is excluded from 
the ground state — this enables us to obtain a reliable 
magnetic phase diagram within finite-size studies by sim- 
ply looking at Stot without worrying about the existence 
of Stot = ferromagnetic-like states. We stress here that 
this procedure is allowed because the SP state may be 
regarded as ferromagnetic: if the SP state were a dis- 
tinct state, then the exclusion of SP states would lead to 
a missed phase transition. 

Here we illustrate this with an appropriate boundary 
condition that makes the F state take a closed-shell con- 
figuration. We present the phase diagram thus obtained 
for the models employed in section III, i.e., ID t-t' Hub- 
bard model, Tasaki's model, and 2-leg Hubbard ladder. 
In all cases, the results accurately coincide with those 
obtained by DMRG calculation for much larger systems 
with OBC. In other words, the result converges to the 
thermodynamic limit rapidly when we concentrate on the 
F states. We may expect that this method should be ap- 
plicable to the determination of the phase diagram of 2D 
or higher dimensional systems. 



A. t-t' Hubbard model 



We first observe that, in the one-electron energy band 
of the t-t' Hubbard model. 



e{k) = -2t cos k + 2t' cos(2fc), 
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the band minimum at fc = splits into double minima 
for t' > 0.25 as t' is increased. Hereafter in this section, 
we set t — 1. 

We first look at the case of t' = 0.2 < 0.25. The band 
has a single minimum, so that an even number of fully- 
polarized electrons take a closed-shell configuration for 
all densities n < 1 if we assume APBC to put the k- 
points symmetrically about fc = 0. In Fig. we show 
the exact diagonalization results for 4, 6, 8, or 10 elec- 
trons in a 12-site system (0.33 < n < 0.83), and 4, 6, or 
8 electrons in a 10-site system. We have also plotted the 
DMRG result for the phase boundary obtained by Daul 
and Noack for a 50-site system with OBCLj'til. We can 
immediately see that the results obtained with APBC for 
fairly small systems is close to that obtained for larger 
ones. 

If we move on to the case of a larger t' — 0.8, an even 
number of fully-polarized electrons only take a closed- 
shell configuration for APBC when the Fermi energy 
is higher than ~2t + 2t' . When the Fermi energy is 
lower, the situation depends on how the fc-points are 
located around the double-minimum dispersion in the 
given boundary condition, so we concentrate on the for- 
mer case. In Fig. we show the results for 8 or 10 
electrons in a 12-site system and 8 electrons in a 10-site 
system, we can see that the result for small systems with 
APBC is very close to the DMRG result for a 50-site 
system. 



B. One-dimensional Tasaki's model 

The phase diagram against U and (3 (that controls 
the dispersion of the band) of the ID Tasaki's model Lb. 
chain of triangles) was obtained by Sakamoto and KuboE£l 
for various band fillings {n — 1/2,1/4,3/8) by means 
of the DMRG method with PBC. As they mention, the 
phase boundary between S'tot = •S'max and ^tot = varies 
strongly with the system size for small systems, so that 
a reliable phase diagram can be obtained only by going 
to sufficiently large (L = 32) systems. 

In Fig. we show the phase boundary determined 
by exact diagonalization for n = 1/2 and system size as 
small as 8 or 12 sites with APBC. Here we fix a = 2, so 
that the lower band becomes flat ioi (3 — 2. Again, the 
result excellently agrees with that of 32 sites with PBC. 



C. Two-leg Hubbard ladder model 

As mentioned in Sec. IIIC, KohncS obtained the 
phase diagram of ?7 = oo Hubbard ladder by means of 
the DMRG method. According to his results, there is a 
wide region of partially ferromagnetic phase in the neigh- 
borhood of paramagnetic phase. 

Let us consider whether we can reproduce the phase 
boundary between non-magnetic phase and partially po- 



larized phase by assuming APBC for small systems. We 
consider the case of 6 and 8 electrons on a 6 x 2 lattice 
and 4 and 6 electrons on a 4 x 2 lattice. We performed 
calculation for 0.125 < t±/{t±+t) < 0.909, and found the 
ferromagnctic-non-magnetic transition only for the case 
of 8 electrons on a 6 x 2 lattice, around t±/LU_+t) ~ 0.45. 
This is consistent with Kohno's the resultlla. 

All these results indicate that the phase boundary be- 
tween the ferromagnetic and the antiferromagnetic states 
is reasonably close to the thermodynamic limit for small 
systems when a care is taken (i.e., letting fully-polarized 
electrons take closed shell configurations). 



V. SUMMARY 

To summarize, we have studied the relation between 
the fully-polarized ferromagnetic state and the spiral 
spins state, a spin-singlet spin state that has a spin corre- 
lation length as large as the system size, which accompa- 
nies the fully-polarized ferromagnetic state in a number 
of electron correlation models. As a typical example, we 
calculate the energy of the SP state and the F state for 
one-dimensional U — oo t-t' Hubbard model, It suggests 
that the SP state is also the ground state in the ther- 
modynamic limit when the F state is the ground state. 
Following the argument by Koma and TasakicSl, we have 
also indicated how a symmetry becomes broken in the SP 
state in the thermodynamic limit, where some 'low-lying 
states' become hybridized. 

We have then characterized the SP state, an itinerant 
magnetic state, by calculating the one-particle spectral 
function. The result is interpreted in a picture in which 
holes move almost freely in a twisted spin configuration. 
We have also calculated the dynamical spin and charge 
correlation functions to find that their behaviors in the 
SP state are similar to those of the F state. 

From these we have conjectured that we should regard 
the SP state and the F state are equivalent in the ther- 
modynamic limit, even though the SP state (spin-singlet) 
and the F state (fully-polarized) are opposite extremes 
in terms of Stot- We have then shown that the mag- 
netic phase diagram can be determined accurately from 
finite-size studies by taking appropriate boundary condi- 
tion (that depends on the number of electrons which are 
accommodated in the shells in the non-interacting case) 
that pushes up the energy of the SP state. This enables 
us to concentrate on the F states, i.e., to simply look at 
the change in S'tot • This is permissible since the SP and F 
states are regarded to be equivalent, so that we can con- 
centrate on either of them. We have obtained the phase 
diagram in this way for ID t-t' Hubbard model, Tasaki's 
model and 2-leg Hubbard ladder. We have found that the 
phase boundary between the ferromagnetic state and the 
non-magnetic state determined in such a way accurately 
coincides with those obtained by the DMRG calculation 
for much larger systems. Since the method does not de- 
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pend on the dimensionality, we may expect that it should 
be applicable to the determination of the phase diagram 
in two or higher dimensional systems. 
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FIG. 2. The DMRG result for the energy differ- 
ence between the SP state and the ferromagnetic state, 
Ai? = Esp — Ef), in ID t-t' Hubbard model as a function 
of inverse system size 1/L for the density of electrons n = 0.5 
with the exchange interaction J = (solid lines), J = 0.2 
(dashed lines), or J = —0.2 (dotted lines). 



FIG. 3. The energy difference between the SP state and 
the ferromagnetic state, AE = Esp — Ep), as a function of 
inverse system size 1/L for the density of electrons n — 0.5 
with the ordinary OBC (dashed line) or the boundary condi- 
tion that omits the hopping at the both ends of the system 
(solid line). 



FIG. 4. The spin-spin correlation function for 10 electrons 
on a 12-site t-t' Hubbard model with (7 = 40 and the periodic 
boundary condition. 
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FIG. 5. Single particle spectral function for 10 electrons 
on a 12-site t-t' Hubbard model with U = 40. Solid lines 
denote the electron-removal spectrum, while dotted lines the 
electron-addition spectrum. Dashed curves indicate the en- 
ergy dispersion defined by eq. ([l^). 



FIG. 6. A plot similar to Fig. |^for 6 electrons on a 12-site 
t-t' model with U = 4 (solid line) and 6 (dashed line) . 



FIG. 7. A similar plot to Fig. |^for 6 electrons on a 12-site 
t-t' model with U — 6. 



FIG. 8. The ground state energy as a function of 1/U for 
6 electrons on a 12-site t-t' model. 



FIG. 9. The dynamical spin(a) and charge(b) correlation 
functions for the SP state (dashed line) and the F state (solid 
line) for 10 electrons on a 12-site ID t-t' Hubbard model with 
U = 40. 



FIG. 10. The spin-spin correlation function for 10 elec- 
trons on a 12-site ID Tasaki's model for U = 4, 6. 



FIG. 11. A plot similar to Fig. |^ for 6 electrons on a 
12-site ID Tasaki's model for U — 12. We have superposed 
the results when we take 7^ = cl- and 7^ = cl , . 
Dashed curves indicate the dispersion defined by eq. (|l2[). 



FIG. 1. The DMRG result for the spin-spin correlation 
function for 24 electrons on a 32-site t-t' model with OBC. 



FIG. 12. The spin-spin correlation function for 10 elec- 
trons ona8x2t/ = 00 Hubbard ladder. 
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FIG. 13. A similar plot to Fig. |5| for 14 electrons on a 
8x2 ladder. Dashed curves indicate the dispersion defined 
byeq. 

FIG. 14. Phase diagram of ID t-t' Hubbard model for 
t' = 0.2 obtained by exact diagonalization for a 12-site sys- 
tem(diamonds) pt|lO-site system(circles) . A DMRG result of 
Daul and NoackEj'tJ is superposed (a full line). 



FIG. 15. A plot similar to Fig. |l|, for the t-t' Hubbard 
model with t' = 0.8. 



FIG. 16. Phase diagram against U and [3 for ID Tasaki's 
model with n — 1/2 for a 12-site system(diamonds) or 8-site 
systep^circles) with PBC. A DMRG result by Sakamoto and 
Kubdij are superposed (a solid line) . 
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